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ABSTRACT 

The Hamburg/ESO quasar HE 1113— 0641 is found to be a quadruple gravitational lens, based 
on observations with the twin 6.5m Magellan telescopes at the Las Campanas Observatory, and 
subsequently with the Hubble Space Telescope. The z$ — 1.235 quasar appears in a cross configuration, 
with i' band magnitudes ranging from 18.0 to 18.8. With a maximum image separation of 0'.'67, this 
is the smallest-separation quadruple ever identified using a ground-based optical telescope. PSF 
subtraction reveals a faint lensing galaxy. A simple lens model succeeds in predicting the observed 
positions of the components, but fails to match their observed flux ratios by up to a magnitude. We 
estimate the redshift of the lensing galaxy to be zl ~ 0.7. Time delay estimates are on the order of a 
day, suggesting that the flux ratio anomalies are not due to variability of the quasar, but may result 
from substructure or microlensing in the lens galaxy. 

Subject headings: gravitational lensing — quasars: individual(HE 1113— 0641) 



1. INTRODUCTION 

A quasar lensed by a galaxy along the line of sight 
presents opportunities to study both the gravitational 
potential that is doing the lensing and the quasar that is 
lensed. In particular the inevitable microlensing of the 
quasar by stars in the lensing galaxies permits estimates 
of the stellar (as oppos ed to dark) matter surface density 



of the lensing galaxy (Schechter & Wambsg anss 20041 
and of the size of the quasar accretion disk ([Pooley et a' 



2007 Poindextcr ct al. 



( |Keeton et al.||2006p . 



20071 and emission line region 
Quadruply imaged systems are 
especially useful for such purposes, giving multiple lines 
of sight through the lens and additional constraints on 
the lensing potential. 

But the rate at which new systems are discovered has 
leveled off to something like a half dozen per year, with 
the most recent discoverie s coming largely fro m the Sloan 
Digital Sky Survey (e.g. |Oguri et aL|2007[ ). Their se- 
lection algorithm is incomplete at separations below one 
arcse cond and worse for quadruple systems than for dou- 
bles flOguri et aLl|2006|j . Of the new SDSS lenses, just 
three are quadruple. No comparable survey has been 
carried out in the southern hemisphere. 

We report here the discovery of a new quadruply im- 
aged system, HE 1113— 0641. It was observed as part of 
a program to identify len sed systems among th e Ham- 
burg/ESO quasar survey ( Wisotzki et al.|[2000[ ). While 
included within the SDSS survey area, it has not, until 
the present work, been identified as a lensed system. 
In f|2j we report the observations made using Magellan 

1 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with program #9744. 

2 This paper includes data gathered with the 6.5 meter Magellan 
Telescopes located at Las Campanas Observatory, Chile. 

3 Massachusetts Institute of Technology, Department of Physics 
and Kavli Institute for Astrophysics and Space Research, 70 Vassar 
St Cambridge, MA 02139; jeffb@space.mit.edu 

Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 
14482 Potsdam, Germany 



and the Hubble Space Telescope, ^describes our analysis 
of the data. In Q we describe a simple model of the 
lens, and in Sj5] we make a rough estimate of the lens 
redshift. In <f6Twe discuss the conclusions we can come 
to regarding HE 1113— 0641. Throughout this paper, we 
adopt a "concordance" cosmology with fl M = 0.3, = 
0.7, and h = 0.7. 

2. OBSERVATIONS 

HE 1113— 0641 was originally discovered to be a z — 
1.235 quasar in the Hamburg/ES O digital objective 



prism survey (Wisotzki et al. 2000). Based on its red- 



shift and apparent magnitude B = 17.01, it was found 
to have a relatively high lensing probability, and was se- 
lected for a follow-up observation. 

We observed the object in early 2002 and early 2003 
in the Sloan g' and i' bands using the Magellan 6.5 me- 
ter telescopes. In Autumn 2003, we observed it in three 
bands using the ACS and NICMOS detectors on the Hub- 
ble Space Telescope. These observations are tabulated in 
Table EJ 

2.1. Magellan 6.5 Meter Imaging 

HE 1113—0641 was identified as a quadruple gravita- 
tional lens on 2002 February 16 using the Baade 6.5 m 
telescope at the Las Campanas Observatory. Seven 60- 
second exposures in each of the Sloan i' and g' bands, 
and a single exposure in the r' band, were taken using the 
Raymond and Beverly Sacklcr Magellan Instant Camera 
(MagIC), a 2048x2048 pixel camera with a 2.4 arcminute 
field of view. The average seeing varied from ff.'43 in i' 
to 0'.'50 in r' to 0'.'52 in g' . Because of the combination 
of mediocre seeing with the existence of only a single im- 
age in the r' band, and the absence of any corresponding 
images in the 2003 dataset, we did not carry out any 
analysis in r' . 

Second-epoch images were obtained on 2003 January 
26, again using MagIC, which had meanwhile been 
moved to the Clay telescope, 60 meters to the north- 
west. The three i' band images had an average seeing 
of 0'.'33, while the single g' band image had a seeing of 
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TABLE 1 
Observations of HE 1113- 



0641 



Date 



Instrument Filter Exposure Seeing 



2002 February 16 


MagIC 


g' 


7 


X 


60 sec 


Of.' 52 






i' 


7 


X 


60 sec 


0'.'43 


2003 January 26 


MagIC 


g' 


1 


X 


60 sec 


0'.'47 






i' 


3 


X 


120 sec 


0'.'33 


2003 November 06 


NICMOS 


H 


3 


X 


640 sec 










1 


X 


704 sec 




2003 November 07 


ACS/WFC 


V 


3 


X 


120 sec 










2 


X 


480 sec 








1 


3 


X 


85 sec 










2 


X 


346 sec 










1 


X 


370 sec 





0"47. 

The data were bias-corrected, flattened, and combined 
using standard techniques. The stacked 2003 i' band 
image may be seen in Figure [l] 

2.2. Hubble Space Telescope Imaging 

On 2003 November 6 and 7, HE 1113-0641 was ob- 
served using both the NIC2 camera of the Near Infrared 
Camera and Multi Object Spectrometer (NICMOS) and 
the Wide Field channel of the Advanced Camera for Sur- 
veys (ACS). The NICMOS images had 256x256 pixels 
and a 19'.'2 field of view, while those from the ACS were 
significantly larger, with 4096x4096 pixels filling a 3.4 
arcminute field of view. Three filters were used, F160W 
with NICMOS and F555W and F814W with the ACS 
(hereafter H, V, and /, respectively). Because of the 
diffraction-limited quality of the images, they were not 
well sampled, with the width of PSF ranging from 1.5 
pixels in H to 2.1 pixels in /. 

We u sed the Multidrizzle program of [Kockcmoc r et aT] 
(2002b I, version 2.2, to register the ACS images, clean 
them of cosmic rays, and drizzle them onto a single im- 
age per filter. The drizzling process also corrects for ge- 
ometric distortion arising from the design of the camera. 
We likewise drizzled the NICMOS im ages into a single 
i mage u sing the procedure detailed in Kockcr noer et al.] 
fl2002a| ) 5 . 

The drizzled ACS and NICMOS images of 



HE 1113—0641 may be seen in Figure [2] 
3. ANALYSIS 

The small separation of this lens, combined with the 
relative faintness of the lensing galaxy, complicated the 
task of disentangling the four quasar components and 
the galaxy, particularly for the ground-based data. To 
address this issue, we fitted for the relative positions and 
brightnesses of the objects using an iterative process. 
First we performed a fit to each image for the relative 
positions of the objects, then averaged the values and 
repeated the fit, holding constant the relative positions, 
to determine the photometry. 

W e used a variant of the DoPHOT photometry pack- 
age ( Schechter et al.||l993 l called Clumpfit to carry out 
the fits, using empirical PSFs provided by field stars for 
the quasar components and a pseudo-gaussian profile for 
the lensing galaxy. (Though this is not a physical profile 



5 The HST Dither Handbook jKoekemoer et al.||2002a| is 
available at: 

http:/ /www. stsci.edu/hst/wfpc2/documents/dither _handbook.html 



choice, we found that the choice of galaxy profile had a 
negligible effect on the goodness of the fit.) We also used 
DoPHOT to obtain astrometry and aperture photome- 
try for several other stars in the wider-field (ACS and 
Magellan) images. 

3.1. Magellan data 

A fit consisting of four empirical PSFs (provided by 
a field star) was made to the stacked i' band Magel- 
lan images. This came to a total of 13 free parame- 
ters: two-dimensional position and normalization for all 
four objects, and the sky level. It was clear from the 
residual images that a small amount of leftover flux re- 
mained; indeed, that the four point sources had been 
over-subtracted in an attempt to compensate (see Figure 
[3| . So a second fit was attempted using a model with four 
point sources and a circularly symmetric galaxy; how- 
ever, there was not enough galaxy flux for the fit to dis- 
tinguish between radial profiles or converge upon a scale 
size. We therefore chose a fixed- width pseudo-gaussian 
profile for the galaxy and repeated the fit, which now had 
16 free parameters. The width of the g alaxy w as set to 
the width chosen for the ACS data (see £ 3.2.1 1, suitably 



broadened by the PSF. This fit was able to account for 
the leftover flux. 

The relative astrometry resulting from the 16- 
parameter fits in i' and other filters is listed in Table 
[2| The relative positions of the four quasar components 
were then weighted by the inverse of their uncertainties 
and averaged, yielding final values, which may be found 
m Table [5] In the case of the len sing g alaxy, only the 
HST positions were averaged (see S 3.2. 1 1 . 

Once the relative astrometry had thus been deter- 
mined, we repeated the fit with fixed relative positions 
and a fixed galaxy width. The results of this fit may be 
seen in Table [4] and residual images are in Figure [2] 

For the g' band images, which had poorer seeing, the 
quasar components were too blurred for a successful fit 
until relative positions were fixed. The residual images 
of these fits may also be seen in Figure [2] There was no 
indication in the residuals of leftover flux indicative of a 
lens galaxy, so we conclude that we have not detected it 
in the g' band. 

Aperture photometry was also obtained for several field 
stars, including those used as model PSFs. The positions 
of these stars, labeled a through e in Figure [T] may be 
found in Table |3l and their magnitudes are listed in Table 

® 

To enable absolute flux calibration, aperture photom- 
et ry was obtained for standard stars from the sample 



of |Smith et~aLl ( |2002[ ). For the 2002 data we used PG 



1047+003A, and forthat of 2003 we used RU 152. We ap- 
plied a first-order correction for atmospheric extinction 
when calculating the z eropoints, using ext inction coeffi- 



cients from Table 4 of |Smith et al.| ( 20071 1 . It is worth 
noting that the u'g'r'i z' system is a broadband approx- 
imation to the (monochromatic) AB magnitude system, 
and is given by 

m = -2.5 log ^.^ff -48.60, (1) 

which is calibrated by synthetic spectra of BD+17 4708. 
The u'g'r 'i' z' deviates from the true AB system by less 
than 5% ( |Smith et al.|2002[ ), and is presented in Table [5] 
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TABLE 2 

Relative Astrometry of HE1113-0641 a 





B 

Aa cos (5) 


AS 


C 

Aa cos (8) 


AS 


D 

Aa cos (<5) 


AS 


G 

Aa cos (<5) 


A<5 


i' (2002) 
i' (2003) 

V 

I 

H 


-C/515 
-0'.'517 
-0'.'518 
-C/519 
-0'.'518 


+0'/428 
+0'.'424 
+0'/424 
+0'/422 
+0'/425 


-0"515 
-0'.'523 
-0'.'523 
-0'/523 
-0'/522 


-0'.'091 
-C/086 
-0'.'085 
-C/083 
-C/083 


-0'.'148 
-0'.'149 
-0'.'152 
-0'.'152 
-0'.'146 


+0'.'433 
+0'/432 
+0'/427 
+0'/429 
+0V429 


-0'.'431 
-0'.'422 

-0'.'320 
-0'.'308 


+0"188 
+0^134 

+0"145 
+0"169 



a All positions are given relative to component A. 



e ■ 

- 




Fig. 1. 
on a side. 



. — 2003 i' band image of the HE 1113— 0641 field, trimmed slightly and binned to 0'.'276 per pixel. The image is 2.2 arcminutcs 
i. The quasar is circled, and the five field stars a through e are labeled. North is up; east is to the left. 
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TABLE 3 
Field Stars a through E a 





Aa cos (<5) 


AS 


a 


+40'.' 04 


-6&/50 


b 


-5'.'70 


-79'.'16 


c 


+66'.' 47 


-29'.'91 


d 


+85'.' 97 


-29'.' 07 


e 


+44'.' 84 


+39'.' 90 



a All positions are given relative 
to component A. 



We estimate the uncertainty in the relative photome- 
try to be 0.1 magnitudes in g' and 0.05 magnitudes in 
i' . Absolute photometry is less certain, with error bars a 
factor of 1.4 larger. With these uncertainties the data are 
consistent with a slight overall brightening of all four im- 
ages between 2002 and 2003, probably caused by intrinsic 
variability of the quasar. However, they fail to convinc- 
ingly demonstrate uncorrelated changes in the flux ratios 
over time, even when combined with HST data; such 
variations might have been indicative of microlensing. 

Despite our use of fixed positions for the g' band im- 
ages, there are inconsistencies in the <?' — i' colors of the 
quasar components. It is likely that these are due to the 
confusion caused by mediocre seeing in the g' band. 

There can be little doubt that the lensing galaxy has 
been detected in the i' band in both data sets. However, 
its size and shape remain poorly constrained. By fitting 
a fixed pseudo-gaussian to both epochs of data, we were 
able to estimate its i' band flux, but with substantial un- 
certainty (0.2 magnitudes of difference between epochs). 
We were able to determine the position of the galaxy us- 
ing HST data, but its size and shape remained elusive 
(see jjO). 



3.2. Hubble Space Telescope data 

The Hubble Space Telescope images did not suffer from 
inadequate seeing, but rather from undersampling of the 
PSF, leading to complications in the interpolation of em- 
pirical PSFs. We therefore resampled the ACS images to 
a scale of 0'.'03 per pixel, and the NICMOS image to a 
scale of 0'.'0375 per pixel, when drizzling. 

3.2.1. ACS 

Since the ACS PSF is known to vary across the field of 
view and also with time, we searched the HST archive 
for images with a suitable PSF star located close to the 
position of HE 1113— 0641 on the chip, and obtained at 
a time close to 2003 November 7. In the V band, we 
used a field obtained on 2003 October 7 6 , and in I we 
used a field obtained on 2003 November 25 7 . In order to 
minimize differences in the PSF caused by the drizzling 
process, we used the same Multidrizzle process on these 
images as on the HE 1113— 0641 images. We chose PSF 
stars that were close to the correct position on the chip, 
and were not saturated. In both cases, these stars were 
about 1 magnitude fainter than image A. 

The fits proceeded as they had in the case of the Mag- 
ellan images. There was appreciable leftover flux in the I 

6 The exposure was associated with program #9756, and started 
at 4:08 AM. 

7 The exposure was associated with program #9822, and started 
at 10:07 PM. 



band, concentrated in the center of the lens system (see 
Figure [3]) . Since the noisy residuals of the quasar com- 
ponents again prevented a measurement of the lensing 
galaxy's radial profile or scale size, we fit it as a circular 
pseudo-gaussian with a fixed width of 0'.'35 (broadened 
slightly by the ACS PSF). The width was chosen by in- 
spection of the residual image, since the choice had little 
to no effect on the goodness of fit parameter. No sign of 
the galaxy was visible in the V band residual image. 

Both of these fits were repeated once we had deter- 
mined and fixed the relative positions. The resultant 
residual images are visible in Figure [2] The magnitudes 
were calibrated using AB zeropoint keywords from the 
HST data headers. The HST broadband flux calibration 
is bas ed on syntheti c spectra of four primary white dwarf 
st ars (Bohlin 1995), and agrees with the AB zeropoints 
of jSmitiTer a". | ( [2a02| to within 3% QBohlin fc Gilliland 
2004[ ). The photometric data are presented in TablcEJ 

Astrometric measurements were made using DoPHOT 
on the / band image. A plate solution was found using 
sixteen USNO-B stars. This solution gives the position 
of component A as (ll h 16 m 23*56, -06° 57' 38"6; J2000) 
to a precision of 0*01 in a and O'.'l in S. 

3.2.2. NICMOS 

The analysis of the NICMOS data was similar to that 
of the ACS data. A 2003 September 5 observation of SA 
107-626 provided a model PSF 8 . In this case, the PSF 
star was much brighter than the quasar components. 

Leftover flux from the lensing galaxy was visible in the 
H band image, as it had been in the / band image (see 
Figure [3J. We again modeled it as a pseudo-gaussian, 
with the same fixed width (broadened slightly by the 
NICMOS PSF). 

A second fit was performed with the relative positions 
of the quasar components fixed to the averaged value. 
The photometry that resulted from this second fit was 
calibrated to the AB system using keywords from the 
data headers, and is visible in Table [3] The residual 
image may be seen in Figure [2] 

4. MODELING THE LENS 

Using the Lensmodel software of Keeton ( |2001[ ), we 
modeled the lensing galaxy as a singular isothermal 
sphere model plus external shear. This model consists 
of a projected 2-dimensional lensing potential given by 



* 2D (0) =br- Vcos2(</>-0 7 ), 



(2) 



where b is the monopole Einstein radius in arcseconds, r 
and <fi are the radial and angular components on the sky 
of the position vector 6, and 7 and <^> 7 are the strength 
and direction of the external shear. Note that in this con- 
vention the position of a companion mass causing a shear 
would be along the ±</> 7 direction. (No such perturber 
is observed in this case, consistent with the apparently 
small shear strength.) This model has seven free param- 
eters, and was constrained by the averaged positions for 
the four components and the galaxy. 

We found the best model to have a radius of 0'.'332, 
with 7 = 0.04 and </> 7 = 37.7 degrees east of north. 

8 The exposure was associated with program #9875, and started 
at 4:56 AM. 
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TABLE 4 

-0641 PHOTOMETRY a 





g> (2002) 


9' (2003) 


i' (2002) 


i' (2003) 


V 


1 


H g> 


- i' (2003) 


V — I 


I - H 


A 


18.37 


18.19 


18.02 


17.96 


18.33 


18.32 


18.25 


+0.23 


+0.00 


+0.07 


13 


18.28 


18.24 


18.09 


18.02 


18.40 


18.35 


18.27 


+0.22 


+0.05 


+0.08 


C 


18.53 


18.39 


18.46 


18.37 


18.64 


18.61 


18.74 


+0.02 


+0.03 


-0.13 


D 


18.91 


18.91 


18.79 


18.76 


19.06 


19.01 


18.92 


+0.15 


+0.06 


+0.08 


G 






22.36 


22.17 




22.47 


21.05 






+1.42 


a b 


20.63 


20.70 


18.08 


18.05 


20.16 


17.86 




+2.65 


+2.30 




b 








18.39 


20.66 


18.04 






+2.62 




c 


19.29 


19.34 


16.62 


16.62 


18.91 


16.39 




+2.72 


+2.52 




d 


18.78 


18.77 


18.36 


18.36 


18.69 


18.45 




+0.41 


+0.24 




c 


22.51 


22.59 


20.01 


20.02 


22.06 


19.81 




+2.57 


+2.25 





All magnitudes are in the AB photometric system (see \ 3.1 1. b a through e are field stars 




Fig. 2. — Top row: Stacked images of HE 1113-0641 from the HST(V, I, and H) and Magellan (g' and i'). Magellan data are from 
2003. The images are displayed with logarithmic stretch. Bottom row: Residual images after subtraction of the best model. The images 
are in a linear stretch from — 20cr to 20<r, where a is the sky noise. All images are 4''0 on a side. 



The source position was predicted to be (Aacos(<5), 
A<))=(-0'.'308, 0'.'151) relative to the position of com- 
ponent A. Table [5] contains a summary of the model's 
predictions, compared to observed data. 

The model fits the component positions very well, even 
with the tight constraints provided by the HST. It does 
not, however, correctly predict the flux ratios. This is 
in keeping with experience; opt ical flux anomalies can 
be seen in many lensed quasars ( |Keeton et aT|| 2003). In 
particular, the D component, a saddle point image, is 
predicted to be brightest, but is observed to the be the 
faintest, too faint by as much as a factor of 2.5. The 
predicted position of the lensing galaxy is O'.'Ol from the 
observed position. This is well within the estimated mea- 
surement error. 

Finally, the model allows us to predict the time de- 
lays between the images, given a lens redshift. We used 
zl = 0.7, as estimated in <35| This yields the predicted 
time delays seen in Table [5T with the maximum delay 
being ~1.5 days. We also calculated predicted time de- 
lays for z L = 0.4 and z L = 1.0; these changes reduced 
and increased (respectively) the time delays by a factor 
of ~3. This strong dependence suggests that a measure- 
ment of the time delays might constrain the lens redshift; 
however, the unknown radial mass profile of the lensing 
galaxy is likely to have a similarly strong effect on time 



delays. 

By way of comparison, HE 0435— 1223, which has a 
shape similar to that of HE 1113—0641 but a larger im- 
age separation, has a m aximum time delay of two weeks 
( |Kochanek et al.||2006[ ). 



5. ESTIMATING THE LENS REDSHIFT 

In order to estimate the redshift of the lensing galaxy 
and determine what its I — H color could tell us about its 
morphology, we combined the results of our lens model 
with properties of typical galaxies. 

From the lens strength b = 0'.'332 we found the line-of- 
sight velocity dispersion of the lens using 



D LS 4ttct 2 
D s c 2 



(3) 



(Narayan & Bartelmann 1996 ), where D^s is the angular 
diameter distance from the lens to the source, and Ds 
is the angular diameter distance from the observer to 
the source. These distances depend on both the source 
redshift zg — 1.235 and the unknown lens redshift zl. 

By combinin g this equation with the Faber-Jackson 
relation ( |Faber fc Jackson|1976|) for e lliptical galaxies, or 
the Tully-Fisher relation ( fluffy &z Fisher|1977| ) for spiral 
galaxies, we generated a predicted observed magnitude 
for each filter as a function of lens redshift. The Faber- 
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TABLE 5 

ASTROMETRY k, LENS MODEL 





Observed a 


Predicted 


Observed b 


Predicted 


Predicted 




Ao cos (<5) 


A<5 


Act cos (<5) 


A<5 






Time Delays c 


A 


0" 


0" 


-0'.'002 


+0'.'001 


+16.0 


+12.6 





B 


-0'.'518 


+0"424 


-0'.'518 


+0'.'423 


+15.7 


+15.7 


7.9 


C 


-0'.'522 


-0"085 


-0'.'524 


-0'.'084 


-10.2 


-9.6 


36.9 


D 


-0'.'150 


+0^'429 


-0'.'150 


+0'.'428 


-8.6 


-16.7 


16.4 


G 


-0'.'313 


+0"158 


-0'.'314 


+0'.'148 









a Weighted average of positions from i' band data and HST data. b Flux ratios from H 
band data; normalized so that component B matches the model. Negative magnification 
denotes saddle point images. c In hours. 





Fig. 3. — Residual images of HE 1113— 0641 in the redder filters after only four point sources, and no central galaxy, have been fit and 
subtracted. Leftover flux from the lensing galaxy may be seen near the center. Stretch and size are equal to those of the second row of 
Figure [2] 

Jackson relation is given by 

M T (B) = -19.4 + 5 log /i- 10(logcr-2.3) 
( de Vaucouleurs fc 01son|[l982 ), which becomes 

M B = -18.9 + 51og/i-10(logo--2.3), (4) 



following the approach of Kochanck (19921. We found 
that the median redshift was 0.66, with a 68% confidence 
interval of [0.41,0.88]. This is consistent with the results 
of the Faber-Jackson method. 



after applying the extinction correction Bt 
( de Vaucouleurs et al.|[l976 l, and with Bj 
(IPeterson et al.||1986J. The Tully-Fisher relation is the 
1; 



B T = 0.22 
B T + 0.29 



same, but with a replaced by circular velocity, which for 
an isothermal sphere is just \[2a. 
The predicted magnitude in the i' band is given by 



= M B + DM(z L ) + K Bti ,(z L ) 



(5) 



where DM(zl) is the cosmological distance modulus, 
and KB.i'{ z h) is the generalized K-correction between 
the lensing galaxy's rest-frame B band ma gnitude and 



the o bserved i! band magnitude (see, e.g., |Hogg et al. 



2002J. To calculate the K-correction for an elliptical 
galaxy at each p otential lens redshi f t, we used an SED 
generated by the Bruzual & Chariot ( 2003 ) spectral evo- 
lution code. Our model consisted of a solar-metallicity, 
instantaneous starburst at a redshift of 3.0, followed by 
passive evolution. For a spiral g alaxy, we used an em - 
pirical Scd galaxy spectrum from Coleman et al. (19801, 
redshifted appropriately. 

The observed magnitudes of the lensing galaxy 
matches those predicted for an elliptical galaxy for a 
range of redshifts 0.4 < z < 1.0. A spiral galaxy model 
also matches the H band observations, but would be 
brighter than observed at all redshifts in / by > 1 mag- 
nitude. The galaxy's brightness and colors seem to be 
more consistent with an elliptical galaxy than a spiral. 

In addition to this method, we estimated the proba- 
bility distribution of the lensing galaxy's redshift by cal- 
culating lensing optical depth as a function of redshift, 



6. CONCLUSIONS 

The z s = 1.235 quasar HE 1113-0641 is lensed into a 
cross configuration, with four components ranging from 
18.0 to 18.8 magnitudes in i' . The maximum image sepa- 
ration is 0'.'67. A combination of ground-based and HST 
imaging has yielded reliable astrometry and photometry 
of the four quasar components of HE 1113—0641, as well 
as a good estimate of the position of the lensing galaxy. 
However, we were unable to measure the size or mor- 
phology of the galaxy. By assuming a circular pseudo- 
gaussian shape and fixing a width for the galaxy, we were 
able to estimate its flux in the redder bands. 

A singular isothermal sphere (SIS) model succeeded in 
matching the positions of the quasar components and of 
the lensing galaxy, but was unable to match the observed 
flux ratios. Based on what has been seen with other 
gravitational lenses, it seems likely that this is due to 
perturbations from stellar microlensing or dark matter 
substructure in the lens galaxy. 

The redshift of the lens galaxy remains elusive, but we 
estimate that z L ~ 0.7 ± 0.3, based on its observed flux 
and colors. At this redshift, its velocity dispersion, as 
measured by its lensing potential, would be 180 km/s. 

Because of its small separation, HE 1113— 0641 will 
likely prove difficult to monitor using ground-based tele- 
scopes, and hope for measuring a spectroscopic redshift 
of the lensing galaxy is faint. Nevertheless, it is an inter- 
esting example of a small-separation lens, and may prove 
useful for studies that can take advantage of telescopes 
with very good seeing. 
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